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ABSTRACT 

We present a flux variability study of simultaneous RXTE and EUVE observations of 
the highly variable Seyfert galaxy NGC 4051. We find a strong correlation between 
variability in the EUV and medium energy X-ray bands, indicating that both are sam- 
pling the same power-law continuum. The lag between the two bands is less than 20 ks 
and, depending on model assumptions, may be < 1 ks. We examine the consequences 
of such a small lag in the context of simple Comptonisation models for the production 
of the power-law continuum. A lag of < 1 ks implies that the size of the Comptonising 
region is less than 20 Schwarzschild radii for a black hole of mass > 10 6 Mq. 
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1 INTRODUCTION 

Despite two decades of observations, the mechanism for 
the production of the X-ray continuum of radio-quiet AGN 
remains poorly understood. The power-law nature of the 
continuum implies a Comptonisation origin, possibly up- 
scattering of low energy (optical, UV) photons by a ther- 
mal distribution of energetic electrons (as indicated by the 
OSSE detection of high energy cutoffs in Seyfert spectra 
e.g., Zdziarski et al., 1997). Rapid variability of the contin- 
uum implies that the X-rays are produced in a small region, 
close to the putative central black hole. These conclusions 
are still vague; further understanding of the innermost re- 
gions of AGN requires more detailed knowledge of the size 
and geometry of the X-ray producing region, the source of 
the low energy 'seed' photons and the physical characteris- 
tics of the scattering particles. 

We can begin to address these points by answering two ques- 
tions: 

1. How low (in energy) does the power-law continuum 
extend? 

2. Are there any lags between the variations in the X- 
ray band and lower energy bands? 

Answering the first question will allow us to place an 
upper limit on the energy of the seed photons (and hence 
allow us to constrain the temperature of the source of the 
seed photons). Answering the second question will allow us 
to determine the size and possibly the geometry of the emit- 
ting region. In upscattering of low energy photons by a cloud 
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of Comptonising particles, more scatterings are required to 
produce high energy photons than lower energy photons, 
hence there will be an intrinsic lag between high and low 
energy bands. The size of this lag is determined by the op- 
tical depth and size of the Comptonising region. 

Searches for lags between the X-ray band and lower- 
energy contnuum bands in AGN have proved difficult, due 
to the difficulty of scheduling simultaneous observations be- 
tween space-based and ground-based instruments for the 
long durations required to search for correlated variability. 
In this respect, the best studied radio-quiet AGN to date are 
the Seyfert galaxies NGC 4151, NGC 5548 and NGC 7469. 
The X-ray/UV monitoring of NGC 4151 and NGC 5548 con- 
strained any lags between the UV and X-ray continua to be 
less than 0.15 d (Edelson et al. 1996) and less than 6 d 
(Clavel et al. 1992) respectively. The intense, month-long 
monitoring of NGC 7469 with the Rossi X-ray Timing Ex- 
plorer (RXTE) and IUE (Nandra et al. 1998) revealed a 
complex correlated behaviour. In this case, the UV contin- 
uum led the X-rays by ~ 4 d at the peaks in the lightcurves, 
while there was no apparent lag between both bands during 
the lightcurve minima, suggesting that complex, multi-zone 
Comptonisation models may be required to explain the rela- 
tionship between the X-ray and UV continua in NGC 7469. 
The similarity between X-ray spectral shape and variability 
properties (e.g. the power spectrum, M c Hardy 1988, Edelson 
& Nandra 1998) of Galactic Black Hole Candidates (GB- 
HCs) and AGN suggests that GBHCs may be scaled-down 
analogs of AGN. A lag of ~ 6 ms has been measured between 
medium-energy (5-14 keV) and low-energy (2-5 keV) X-ray 
bands in the GBHC Cygnus X-l (Page 1985), which scales 
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to a lag of order 100 s-10 ks in AGN, if we naively scale 
by the black hole mass (for a range of black hole masses, 

10 6 10 s Mq, as measured by Wandel Peterson & Malkan 

1999). In simple Comptonisation models, we expect lags be- 
tween bands to be energy dependent, with greater lags be- 
ing seen between bands which are more separated in energy. 
This expected dependence has been seen in GBHCs (e.g. 
Cygnus X-l, Nowak et al. 1999), so that if the same ap- 
plies to AGN, we might expect lags between the X-ray and 
EUV bands to be less than those between the X-ray and UV 
bands (i.e. less than 1 day) and more easily measurable by 
a relatively short monitoring campaign. 

NGC 4051 is a nearby (z = 0.0023), low luminosity 
(Lx-ray < 10 42 ergs s _1 ) Narrow Line Seyfert 1 galaxy, 
which display s strong X-ray variability in b oth medium en- 
ergy X-rays (Papadakis & Lawrence 1995) and low energy 
X-rays (Lawrence et al. 1987, M c Hardy et al. 1995) with a 
doubling timescale of less than 1 ks. This strong variability 
makes NGC 4051 an ideal candidate for searching for lags 
between the X-ray band and lower-energy bands. We there- 
fore obtained simultaneous observations of NGC 4051 with 
RXTE and the Extreme Ultraviolet Explorer (EUVE) satel- 
lite in May and December 1996. In this paper we present the 
results of these observations as they relate to flux variability. 
Detailed spectral analysis of the RXTE observations is left 
to a later work. 

In the next section we discuss the observations and 
data reduction method, while in section 3 we present the 
lightcurves and their cross-correlation function. In section 4 
we further constrain any lag by assuming a simple spectral 
model. Finally, we discuss the implications of our results 
for the size of the scattering region in NGC 4051 and for 
Comptonisation models of the X-ray continuum. 



2 OBSERVATIONS AND DATA REDUCTION 

From May 20 to May 28 1996 RXTE observations of typ- 
ical duration ~ 1 ks were carried out approximately every 
half-day and continuous (subject to source visibility) EUVE 
observations were also carried out at this time. Additional 
simultaneous observations occurred on May 6 1996. From 
December 13 to December 16, a long RXTE observation was 
carried out, which overlapped with a continuous three-day 
EUVE observation for ~ 120 ks. 

The RXTE satellite observed NGC 4051 with the Pro- 
portional Counter Array (PCA) and the High Energy X-ray 
Timing Experiment (HEXTE) instruments. The PCA con- 
sists of 5 Xenon Proportional Counter Units (PCUs), sensi- 
tive to X-ray energies from 2-60 keV. The HEXTE covers a 
range of between 20-200 keV, but due to the faint nature of 
the source we only consider the PCA data in this work. 

Discharge problems mean that of the 5 PCUs in the 
PCA, PCUs 3 and 4 are often switched off, so we include 
data from PCUs 0, 1 and 2 only. We extract data from 
the top layer of the PCA using the standard ftools 4.1 
package, excluding data obtained within and up to 20 min- 
utes after SAA maximum and data obtained with earth el- 
evation < 10°. We generate background data for the PCA 
with PCABACKEST v2.0c using the new L7 model for faint 
sources. 

NGC 4051 was observed with the Deep Survey Spec- 



trometer (DS/S) on board the EUVE satellite. The DS/S 
(e.g. Welsh et al. 1990) is equipped with a broad band imag- 
ing detector (covering the 66-178 A band in the Lexan/B 
filter) and three spectrometers covering the 'short' (SW: 70- 
190 A), 'medium' (MW: 140-380 A) and 'long' (LW: 280- 
760 A) EUV wavelengths (Abbott et al. 1997). This config- 
uration allows simultaneous imaging and spectroscopy with 
a spatial resolution of ~ 1 arcmin and a spectral resolution 
of A/AA ~ 200 at the short wavelengths. 

NGC 4051 was detected only in the SW spectrometer 
because the interstellar medium (JVh = 1.31 x 10 20 cm~ 2 
Elvis, Lockman & Wilkes 1989) severely attenuates the EUV 
spectrum for wavelengths longer than ~ 100 A. But despite 
the long (~ 203 ks) total effective exposure, the average 
signal-to-noise ratio (SNR) achieved in the spectrum - per 
1 A bin - is very low (~ 1.5cr in the range 80-90 A) and spec- 
tral data will not be discussed further here. For the purposes 
of this work, we only use the lightcurves produced by the DS 
imaging detector in the energy range 124-188 eV (with the 
Lexan/B filter). 

We extracted the light curves from the DS time-ordered 
event list using the EUVE Guest Observer Center software 
(IRAF/EUV package) and other IRAF timing tasks adapted 
for EUVE data. After correcting the data for instrumental 
deadtime, telemetry saturation, vignetting and eliminating 
intervals of high particle background, the total effective ex- 
posure is 178738 s. We counted the source photons in a circle 
of 2' radius, and we estimated the background in a concen- 
tric annulus with inner and outer radii of 2.7' and 8' respec- 
tively. The extraction region includes more than 98% of the 
DS point spread function (Sirk et al. 1997) and the back- 
ground is generally very uniform in the chosen area. The ef- 
fective area of the DS instrument (25 cm 2 at A = 85 A, Sirk 
et al. 1997), is more than 10 times larger than the spectrom- 
eter effective area at this wavelength and a good detection of 
the source was obtained during each 5544 s EUVE orbit (see 
figures 1 and 2). A complete description of the EUVE ob- 
servation and data analysis can be found in Fruscione et al. 
(1999). Similar lightcurve extractions have been performed 
on other AGNs detected by EUVE (e.g. 3C273, Ramos et 
al. 1997, NGC 5548, Marshall et al. 1997). 



3 TIMING ANALYSIS 

We show the background-subtracted RXTE (2-10 keV) 
and EUVE lightcurves for the May observations in Fig. 1. 
For presentation purposes, we have scaled up the EUVE 
lightcurve by a factor of 100. The bin width is 5544 s, which 
is the orbital period of EUVE. For clarity we do not show 
the May 6th observations (but we note that the 2-10 keV 
and EUVE count rates for this time were ~ 8.0 cts s~ 
and 0.12 cts s _1 respectively). As can be seen from the fig- 
ure, the EUVE count rate varies strongly between less than 
0.05 cts s" 1 and ~ 0.25 cts s" 1 . The 2-10 keV count rate 
also varies greatly between RXTE observations, with a range 
of between 3 cts s _1 and 17.4 cts s _1 . The lightcurves for 
the December observations are shown in Fig. 2, using the 
same bin width and the same scaling factor for the EUVE 
lighcurve including, for completeness, the data obtained by 
both instruments outside the times of overlap. Strong vari- 
ability can be seen, with NGC 4051 reaching a particularly 



© 0000 RAS, MNRAS 000, 000-000 



EUV and X-ray variability of NGC 4051 3 




142 



144 



146 

Days in 1996 



148 



150 



Figure 1. May 1996: 5544 s binned EUVE (scaled by factor 100) and RXTE (3 PCUs, 2-10 keV) lightcurves of NGC 4051 (errors are 
la). 
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Figure 2. December 1996: 5544 s binned EUVE (scaled by factor 100) and RXTE (3 PCUs, 2-10 kcV) lightcurves of NGC 4051 (errors 
are lcr). 
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Figure 3. Z-transformed discrete cross-correlation function of EUVE and RXTE (4-10 keV) lightcurves, combined from separate May 
and December ZDCFs. The ZDCF is binned to 20 ks and errors are la. 



low flux state during the period of overlap. Visual inspec- 
tion of the lightcurves shows a striking correlation between 
the two bands. We can quantify this correlation and search 
for lags between the two bands by carrying out a cross- 
correlation analysis of the lightcurves. In the context of this 
paper we are particularly interested in the relation of the 
power-law component of the X-ray spectrum to the EUV 
emission. Since NGC 4051 is known to display both vari- 
able low energy X-ray abs orption and an iron fl uorescent 
emission line at ~ 6 keV (Guainazzi et al. 1996), we will 
exclude these components and obtain the 'pure' power-law 
contribution to the emission by extracting a lightcurve in 
the 4-10 keV band, excluding 5-7 keV (corresponding to 
PCA channels 12-14 and 21-28 in the current gain epoch 
3). 

We compute the cross-correlation function (CCF) of 
the time series using the Z-transformed Discrete Correla- 
tion Function (ZDCF) method of Alexander (1997), which 
is based on the DCF method of Edelson & Krolik (1988) but 
estimates errors more reliably. We calculate seperate ZDCFs 
for May and December 1996 and bin up the resulting noisy 
ZDCFs (which have different binning) into identical 20-ks- 
wide bins, before combining them by adding together values 
corresponding to the same lag (weighting according to er- 
rors). We show the resulting combined ZDCF in Fig. 3. 

A simple visual inspection of the ZDCF indicates that 
the EUV band leads the X-rays by between and 20 ks. Sim- 
ulations of perfectly correlated lightcurves (with zero lag) 
which use the same sampling pattern as our data show that 
the peak at 20 ks is probably artificial, due to an excess 
of data pairs sampled at that particular lag. Therefore, we 
believe that the most likely lag is within the ks bin. We 
now turn to an alternative technique for constraining the 
lag, which is not affected by the sampling pattern. 



4 A SIMPLE SPECTRAL MODEL 

The simplest explanation for the strong correlation between 
the EUVE and RXTE lightcurves is that both instruments 
are sampling the same continuum, i.e the X-ray power-law 
extends to the EUV band. A detailed spectral analysis of the 
RXTE data will be described in a later work, but here we 
can test the hypothesis that the EUV continuum is an exten- 
sion of the X-ray power-law continuum by fitting a simple 
power-law to the 4-10 keV (excluding 5-7 keV) region of 
the spectrum for the brightest observation in May. We find 
a power law slope (photon index, V) of 2.3±0.1, correspond- 
ing to a 2-10 keV flux of 6.5 x 10" 11 ergs cm" 2 s" 1 . If we 
use these parameters in the pimms v2.3 count rate calcu- 
lator, assuming a galactic absorption of 1.31 x 10 20 cm" 2 , 
we obtain a predicted EUVE DS/S (Lexan/B) count rate of 
0.27 ± 0.1 cts s" 1 . This agrees well with the actual EUVE 
count rate (simultaneous with the 500 s long RXTE obser- 
vation) of 0.24 ±0.05 cts s" 1 . Note that because the EUVE 
count rate can be adequately explained as the extrapola- 
tion of the X-ray power law modified by Galactic absorp- 
tion, there is no requirement for a significant neutral column 
in the AGN host galaxy, a result consistent with previous 
RO SAT and ASCA observations. The lack of any significant 
absorbing column (in addition to Galactic) implies that the 
significant low energy X-ray absorption seen in previous ob- 
servations of NGC 4051 (M c Hardy et al., 1995; Guainazzi et 
al., 1996) must be due to ionised gas (so that hydrogen and 
helium are virtually completely ionised). This is again con- 
sistent with previous observations, which indicate that the 
X-ray absorption is due to an ionised 'warm abso rber'. Pre- 



vious observations of NGC 4051 with EX OSAT flPapadak: 
Sz Lawrence 1995) and recently ASCA ( ]Guainazzi et al 



199f), indicate that the photon index of the X-ray power- 
law is positively correlated with the source luminosity, in 
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Figure 4. Comparison of linear and non-linear models (dotted line and solid line respectively). For clarity the data have been averaged 
(and la errors calculated accordingly) into bins of 0.3 counts s _1 width. See text for details of corresponding fit parameters. 



the sense that the continuum becomes softer as the source 
flux increases. An inspection of the scaled EUVE and RXTE 
lightcurves indicates that this may indeed be the case; there 
is a tendency for the linearly scaled EUVE flux to exceed 
the 2-10 keV count rate when the source is bright, and vice 
versa when the source is dim. If the EUVE count rate can be 
described as a simple function of the RXTE count rate, we 
can fit such a function to the data. The scaled lightcurves 
indicate that the EUVE count rate scales with the RXTE 
count rate in a non-linear way. We suggest a simple function 
of the form: 



Reuve = A _Rp 



, + C 



where -Reuve and Rrxte are the predicted EUVE count 
rate and the RXTE (4-10 keV, excluding 5-7 keV) count 
rate respectively. A, n and C are constants. The error on 
the predicted values is given by: 

n ARrxte 



A Re 



(Re 



C) 



Rrxte 

where AReuve and ARrxte are the errors on the EUVE 
and 4-10 keV count rates respectively. 

For given parameters A, n and C we can calculate the 
predicted value (and error) of -Reuve for each point in the 
RXTE lightcurve and compare with the actual values by 
calculating a \ 2 value, defined as: 

2 (-Reuve. measured — Reuve) 2 

X = x 



/ j \ / >:■ 

" n EUVE,mcasured 



+ AR| UVE 



Where -Reuve, measured and A_R E uve, measured are the mea- 
sured values of the EUVE count rate and associated error 
respectively, and the sum is over the n pairs of data points 



that are measured simultaneously in both bands. By step- 
ping through a range of values of the parameters we can 
attempt to find the set of best-fitting model parameters. 

We use the 4-10 keV band of the RXTE data, (excluding 
5-7 keV), so that we sample only the continuum component 
of the X-ray spectrum. Due to orbital constraints there are 
gaps in the EUVE lightcurve so that not all RXTE data 
points have corresponding simultaneous EUVE data points. 
We therefore bin the lightcurves into identical 1 ks bins (so 
the corresponding time for each bin is the same for both 
the EUVE and 4-10 keV lightcurves). In this way we only 
sample times where the EUVE and RXTE data are simul- 
taneous to within 1 ks. This restriction is necessary because 
large flux changes can occur on timescales of a few ks, so 
the flux in both bands must be measured as simultaneously 
as possible to allow a good comparison between bands. We 
next fit the simultaneous data with our model. 

We first fit a simple linear model to the data (n = 1). 
We obtain best-fit values of A = 0.029=hg;ggt and C = 
— O.006±8:8or, with X 2 = 35.14 for 31 degrees of freedom. 
This model is acceptable, but a non-linear model (allow- 
ing n to be a free parameter) improves the fit {\ 2 /d.o.f. = 
30.00/30), with best-fit values A = 0.012±g;8J|, n = 1.68 ± 
0.55 and C = 0.006±o;jm • Negative and positive values of C 
(as given by the linear and non-linear model fits) correspond 
to constant flux components in the RXTE and EUVE bands 
respectively. Recently NGC 4051 was observed in an ultra- 
low state by BeppoSAX, RXTE and EUVE (Guainazzi et 
al. 1998, Uttley et al. 1999). The BeppoSAX and EUVE ob- 
servations revealed a constant soft component which may be 
associated with the extended emission in the host galaxy im- 
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aged by ROSAT (Singh 1999). This component contributes 
an EUV flux consistent with the constant value estimated 
by our non- linear model fit. Additionally, RXTE and Bep- 
poSAX detected a constant hard component which may 
be due to reflection from a molecular torus. On its own, 
such a component would lead to a negative value of C. 
The exact value of C, which is at present not well con- 
strained, is determined by a trade-off between these two 
constant components. We note also that fluctuations in the 
cosmic X-ray background, which are unaccounted for by the 
PCA background model, could lead to some additional con- 
stant negative or positive offset which will contribute to 
C. However, the good agreement between the simultaneous 
medium-energy X-ray spectra obtained by BeppoSAX and 
RXTE during the ultra-low state (Uttley et al. 1999) shows 
that any offset due to inaccurate background estimation is 
much smaller than any offset due to the torus component. 
Using our data, we cannot formally rule out either the non- 
linear or linear model, however the F-test indicates that the 
non-linear model is better at describing the data at 63% 
confidence. In figure 4 we show a comparison of the linear 
and non-linear models with the data. 

The fact that we can successfully model (reduced x 2 = 
1.0) the relationship between both energy bands to a simul- 
taneous time resolution of 1 ks implies that there is no sig- 
nificant lag between the bands (i.e. any lag is less than 1 ks). 
This result is consistent with the less-than-unity value of the 
zero-lag peak in the lower resolution CCF. We can rule out 
the possibility that our model is so general that it will fit 
any time lag by shifting the lightcurves with respect to each 
other and fitting the model. If we cause the EUVE lightcurve 
to lag the RXTE lightcurve by 1 ks, the best-fitting non- 
linear model yields x 2 /d.o.f. — 65.78/23. If we cause the 
RXTE lightcurve to lag the EUVE lightcurve by 1 ks we 
obtain '/d.o.f. = 39.38/35. If the RXTE lightcurve lags 
the EUVE lightcurve by 2 ks we find X 2 /d.o.f. = 60.78/38. 
Clearly the model does not easily fit other lags, although a 
lag of 1 ks between the X-ray and EUV bands is allowed by 
the model. We might expect the model to fit short lags even 
if the true lag is zero, since both lightcurves are autocorre- 
lated on short timescales (~ks). In any case, it remains true 
that effectively zero lag between the bands is adequate to 
describe the data. 

Recent simultaneous observations of NGC 5548 with 
RXTE and EUVE showed evidence that the 2-20 keV X- 
ray variations lagged changes in the EUV by 35 ks (Chiang 
et al. 1999). Scaling this lag by the black hole mass, using 
the masses of NGC 4051 (1.4 x 10 6 M Q ) and NGC 5548 
(8 x 10 7 Mq ) measured by reverberation mapping (Wandel, 
Peterson & Malkan 1999), we predict an X-ray-EUV lag of 
~ 600 s in NGC 4051, consistent with our upper limit. 



5 DISCUSSION 

We have shown that the variable X-ray and EUV fluxes 
almost certainly come from the same spectral component, 
and th at the best-fitting model to describe the relationship 



served previously by ASCA (Guainazzi et al. 1996 ) . Using 



betweei l the fluxes in both bands is probably non-linear, 
implying that increases in continuum flux are accompanied 
by spectral slope changes. Specifically, the power-law steep- 
ens as the continuum flux increases. This behaviour was ob- 



pimms we can estimate the change in power-law slope as- 
sociated with a doubling of the RXTE count rate assuming 
n = 1.68. The corresponding slope increase is ~ 0.1, a value 
consistent with the results obtained by Guainazzi et al. 

By constraining the lag between the EUV and 4-10 keV 
bands to be less than 20 ks or less than 1 ks (depending on 
whether we take the conservative result of the ZDCF, or 
the more speculative result of the scaling model), we can 
constrain the physical size of the Comptonising region in 
the context of simple upscattering models. In this context, 
the lag between the EUV and 4-10 keV bands corresponds 
to the time taken to upscatter EUV photons to medium 
X-ray energies (regardless of the energy of the initial seed 
photons). We now consider a simple Comptonisation model, 
where the EUV photons originate at the centre of a homo- 
geneous spherical cloud of thermal electrons, and are up- 
scattered on their way through the cloud to produce the 
medium-energy X-ray photons. We use this simple model to 
estimate an upper limit to the size of the X-ray emitting 
region. 

The time spent to upscatter a photon from an energy 
E\ to a higher energy E2 is: 

NX 

iu P = — (1) 

where N is the number of scatterings required to raise the 
energy from E\ to E2 and A is the mean free path of the 
photon between scatterings. The mean free path is given by: 

A = -L- (2) 

where n c is the electron density and ctt is the Thomson 
cross-section (assuming E\ <C m c c 2 ). We can express n c in 
terms of the optical depth of the Comptonising cloud, r and 
the radius of the cloud, R: 

(3) 



(XT R 

Incorporating equations 2 and 3 into equation 1, and rear- 
ranging yields: 

C T tup 



R = 



N 



(4) 



so that substituting our maximum lag for t up sets an upper 
limit on R. The number of collisions required to scatter the 
photons from energy E\ to E2 depends on the electron tem- 
perature, T c (e.g. see Longair, 1992). In the non-relativistic 
regime, the number of scatterings is given by: 

log (E2/E1) 



N = 



log (1 + (4 fcT c /m c c 2 ))' 



(•>) 



The electron temperature is not well constrained, but be- 
cause the power-law continuum extends right across the use- 
ful PCA band (up to 15 keV) with no cut-off at the high- 
est energies, we are justified in assuming a lower limit to 
the electron temperature of ~ 25 keV. OSSE observations 
of Seyfert galaxies show high energy cutoffs in their con- 
tinua which correspond to typical electron tempe ratures of 
~ 100 keV (~ 60 keV in the case of NGC 4151) (ftdziarski 



et al. 1997 ), consistent with our simple non-relativistic as- 



sumption. 

The ratio of photon energies between the two bands is 
> 20, so for an example where k T c = 25 keV, N > 16. For 
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kT c = 100 keV, N > 5. For photon energies E < k T c , the 
optical depth for Compton scattering in a spherical cloud 
can be determined from the spectral slope of the result- 
ing Com ptonised continuum (Posdnyakov, Sobol & Sunyaev 
19830 ' 



0.5 



4 [{a + 3/2) 2 - 9/4] k T c 



where a is the energy spectral index of the Comptonised con- 
tinuum. Using the spectral index of 1.3 measured by RXTE, 
we find t = 1 for k T c = 100 keV, and r = 2.5 for k T c of 
25 keV. Therefore, in the range of temperatures 25-100 keV, 
we expect the ratio t/N < 0.2. Using this limit in equation 
4, we arrive at an upper limit for R assuming t up < 1000 s: 

R < 6 x 10 12 cm. 

A size of 6 x 10 12 cm corresponds to ~ 20 Schwarzschild 
radii for a 10 6 M© black hole. Thus, if the power-law con- 
tinuum is produced by a single central source, the emitting 
region is very close to the central black hole. Alternatively, 
in disk corona models (e.g. see Haardt, Maraschi & Ghis- 
ellini, 1997) the power-law continuum is produced by the 
scattering of thermal photons from an accretion disk in a 
hot corona above the disk. In this case, the continuum may 
originate in many scattering regions spread over the disk, 
so the size constraint we impose may represent the size of a 
typical scattering region. Note that this size limit increases 
by a factor ~ 20 if we instead choose the more conservative 
20 ks upper limit on the lag, inferred from the ZDCF. 

We note here that Monte Carlo simulations of thermal 



Comptonisation in the relativistic regime ( Bkibo et al. 199E ) 
show that if a > 1, r < 0.1 for k T c > 250 keV. In the 
case of very low optical depth to scattering, EUV and X-ray 
photons may undergo only one upscattering from the orig- 
inal seed photons before leaving the Comptonising cloud. 
In this case, upper limits to lags between the two bands 
yield no information regarding the size of the Comptonising 
region. This situation is extreme however, but in the case 
of NGC 4051 we cannot be completely confident in our as- 
sertion of a small Comptonising region until finite lags are 
measured between different energy bands, or a spectral cut- 
off is determined. 

Finally, we note that although our constraint on the 
size of the emitting region in NGC 4051 is model depen- 
dent, simple causal arguments place a model-independent 
upper limit of 3 x 10 13 cm on the separation between the 
EUV and X-ray emitting regions, if the lag between both 
bands is less than 1 ks. Combining this constraint with the 
rapid variability seen in both bands, we can infer that the 
continuum in both bands is emitted co-spatially within a 
region of that size (i.e., 1000 light-seconds). 



6 CONCLUSIONS 

In the introduction to this paper, we raised two key ques- 
tions which would shed light on the nature of the X-ray con- 
tinuum production mechanism in Seyfert Galaxies, namely 
how low in energy does the power-law extend? And are there 
any lags between the X-ray band and lower energy bands? 
We have answered these questions to some extent in the case 
of NGC 4051. We have shown that the EUV variability of 



NGC 4051 can be adequately explained as being due to the 
extension of the X-ray power-law into the EUV band. This 
implies that the source of 'seed' photons required must be 
looked for at lower energies. Done et al. (1990) have shown 
that NGC 4051 does not vary significantly in the optical 
band, so the X-ray power-law must end somewhere in the 
UV, Far UV or in the unobservable, heavily absorbed EUV 
band above the Lyman edge. 

We have shown that there are no significant lags down 
to timescales of 20 ks, and that if the EUV emission is the 
extension of the X-ray power-law (of varying slope or oth- 
erwise), the lag is reduced much further to less than 1 ks. 
This places strong constraints on the size of the Compton- 
ising region, implying that it lies within a few gravitational 
radii of the central black hole for a typical black hole mass 
of between 10 6 and 10 7 Mq , if the X-ray source is central. 

Since NGC 4051 is a low luminosity Seyfert galaxy, we 
might expect that more luminous AGN will contain larger 
continuum emission regions. The lack of strong variability 
in higher luminosity AGN makes the type of study outlined 
in this paper difficult, but the high count rate, high S/N 
observations that will be available with XMM should prove 
excellent for constraining the size of the continuum emission 
region across a range of source luminosities. 
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